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Abstract. Theoverexploitation of groundwater resources can cause several 
problems; one is the use of wells for situations that require a higher de- 
mand on the avail ability of the well reservoir leading to exhaustion. Another 
issue is the lowering of soil compaction with the underlying layers of po- 
rous, due to water extraction capacity in excess of recharge. This phenome- 
non is known as subsidence of the soil. This article shows a general review 
of remote sensing technique, namely DlnSAR - Differential I nterferometry 
Synthetic Aperture Radar applied to studies of land subsidence. From these 
techniques it is possible to seethe lowering of the ground with theoverex- 
ploitation of aquifers. The Dl nSAR technique is independent of other tech- 
niques for measuring subsidence, such as geometric leveling and gravimet- 
ric method, whereas the sensor system Advanced Synthetic Aperture Radar 
orbiter is not located in the soil-building. Other types of phenomena (seis- 
micity, ground works, among other ones) can generate subsidence of land, 
may occur in the same season and region and overlap the overexploitation 
of aquifers. 
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1. Introduction 

The subsidence of the soil is a rheological phenomenon lowering of the 
ground surface in response of changes in the groundwater by reducing the 
level of alti metric ground due to the removal of subsurface support (Cabral 



etal., 2006; Paula, 2005, Cabral etal. 2001 Cabral etal. 2002; Cabral etal. 
2003; Silva, 2004; Chen et al. 2003; Chang et al., 2005; Gonzalez-Moran et 
al. 1999). 

The subsidence of the soil can be triggered by a cause alone or be linked to 
other natural phenomena, namely, theseismicity generated by movements 
of normal faults and transcurrent; earthquakes in subduction zones; in- 
duced earthquakes and volcanic eruptions dams (Takeya et al . 1989); ac- 
commodation sedimentary structures induced by earthquakes collapsing 
caves (karst dissolution); extraction of gas and oil; structural collapse gen- 
erated by engineering works of mine as mine collapsed due to vibrations 
caused by deflagration of explosives and failures structural protection (U.S. 
Army Corps of Engineers, 2002), and overexploitation of water from aqui- 
fers (Cabral, 2001; Santos, 2005). Underground civil engineering (construc- 
tion of underground canals, tunnels, subways, urban drainage tanks) de- 
signed and built incorrectly or without geotechnical surveys, geophysical, 
geological and geo-hidrological appropriate, can generate subsidence. Ac- 
cording to Scott (1979) the phenomenon of subsidence can be seen as a ge- 
ographical expression in the area, the time consumed and human activity or 
natural phenomena isolated inducing subsidence, as shown inTablel 



AREA 

(km 2 ) 


TIME INTERVAL 

(years) 


HUMAN ACTIVITY/ NATURAL PHENOMENON 


10 2 a 10 3 


10 4 


Isostatic adjustments of the crust 


10a 10" 


tectonic plates 


10 a 10" 


Exploitation of water on a large scale 


10"' a 10" 


earthquakes 


10 a 10" 


10a 10" 


Exploitation of oil and gas 


10 u a 10 


10 a 10" 


Extracting solids 


10"' a 10 u 


10a 10" 


Charging the surface, consolidation, hidrocompactation 


10"'a10 u 


10"' a 10" b 


Densification by vibration 


10" a 10" 1 


10" 1 


Compression 



Table L The soil subsidence according area, time, and origin, Scott (1979). 



2. Techniques for measuring subsidence of the soil 

Several independent methods can be found in the literature for measuring 
the seismic subsidence, I ike the gravimetric and magnetic, Laser Scan, aeri- 
al photographs, GNSS (Global Navigation Satellite Systems), which require 
integration with geological information, logs of wells, hydrological maps, 
etc. to allow validation data and isolation of the causes of subsidence 
(Thomas etal. 2004). 



The gravi metric method can be used to measure displacements of the water 
body, subsurface erosion and accommodation (Blackely, 1996; J unior 
Tavares, 2009; J unior Tavares, 2008; Zidarov, 1990), and It can be used in 
cities and rural areas at various scales (1100, 1500, 110000). The gravi- 
metric method not depending on the purchase of images to generate con- 
tour maps showed and a reduced requiring teams with two digital gravi me- 
ters. 

In urban areas the technique for control GNSS altimetry measurement of 
soil subsidence requires simultaneous operation of dual frequency receivers 
same specs screening and calibration of antennas; Either it is necessary 
continuous measurement periods, staff costs, etc., (Leick, 1995; Seeber , 
2003; Romaoetal. 2003; Monico, 2007). 

The GNSS signal interference could be caused by buildings, vehicles vibra- 
tions, and protection of landmarks occupied by receivers hinder systematic 
measurement campaigns (Blachut et al., 1979). 

Alti metric control points need to be outside the area of deformation to ana- 
lyze the stability of the land and appropriate measurement techniques (Mi- 
randa, 2007; Soto, 2006; Burity and Seixas, 2005; Granemann, 2005; Chen 
and Chrzanowski, 1990; Chrzanowski et al. 1985), further defining the geo- 
detic reference system (Sato et al., 2003, Sneed et al. 2002; Sneed et al. 
2003; Dalazoana, 2005, Freitasetal., 1999; Medeiros, 1999). 

Moreover altimetry obtained by the GNSS can be compared to data ob- 
tai ned with the above menti oned techni ques. 

2.1. The PSInSAR technique 

I n the context of interferometric techniques to measure subsidence is im- 
portant to emphasize that the technique PSInSAR represents a paradigm 
for developers academics studyi ng the accuracy of Dl nSAR to measure sub- 
sidence and also create new algorithms capable of detecting PS. 

Researchers I ike Chang and Rizos(2005), Elias et al. (2009), Zhou (2009), 
Froguer et al. (2007) have developed their own mathematical tools and 
methodologies to detect PS from classical DlnSAR repository, obtaining 
good results in the measurement of subsidence on cost/ benefit and accura- 
cy of the order of cm/ year to mm/ year (Duy2009; Thomas etal. 2004). 

Moraetal. (2007) used inthetown of Cardona, inthe central region of Cat- 
alonia, Spain, with DlnSAR, fewer images and precision compatible tech- 
ni que PSInSAR. 

Henriqueset. al. (2011) Describe a study of land subsidence in Lisbon Area 
with Val i dati on Of Psi nSAR Results. 



The technique PSInSAR was created and patented in 1999 by a group of 
electronicsengin eers f r om P ol ytech n i c U n i versi ty of M i I an ( P L I M I ) , I tal y, 
with several applications, among which the measurement of soil subsidence 
(Ferretti, et al. 2001). The group owns the property POLIMI algorithm 
PSInSARTM technology spin-off company in Tele-Rivelamento Europe - 
TER since 2000 on five continents. Mapping subsidence of soil using the 
technique PSInSAR detecting vertical movements can achieve accuracy in 
theorder of magnitude of cm/ year or until mm/ year(Ferreti etal.2001). 

2.2. Principles of the theory of subsidence measurement using 
DlnSAR 

The DlnSAR Remote Sensing technique can be used in the acquisition al- 
timetry and integrating i nterferometry SAR images (Shearer, 1998; Abidin 
et al. 2001 Colesanti et al. 2,005, Matthias, 2006, Matias et al. 2005; Ostir 
and Komac, 2007; Elias et al. 2,009; Chang, 2000; Ferreti et al. 2001 
Lauknes, 2004; Fruneau et al . 2003). 

Figure 1 shows the fundamental DlnSAR geometry. The pointisthe acquisi- 
tion of Om reference Os and the acquisition of repetition. The term B is the 
baseline interferometric with B T component perpendicular and parallel B P . 
The term R M is the slant range in decision Om, Os and Rs in sight. The 
viewing angle 6 is the sensor for one pixel in the soil; d6 is the difference 
between Pi and P2 sight. The sensor is ellipsoidal height h, and the differ- 
ence in altitude between Piand P2 is h P . The angle a to the horizontal sen- 
sor 60 is at zero altitude. 

There are several possible combinations of images for measuring subsid- 
ence in terms of baseline critical, requiring careful analysis to determine a 
correct combination to generate the interferogr am and the degree of coher- 
ence between images. A baseline interferometric too small or too large, can 
disable the generation of the i nterferogram because this factor is based on 
the appropriate acquisition geometry, there is a critical baselineto consider. 
A coherent mix between takes Om and Os determines the phase difference 
pixel by pixel, there are fundamental knowledge and high spatial correla- 
tion ephemeris accurate sensors in M and O s (Lauknes, 2004). 

Differences in the geometry acquisition generate geometric decorrelation. 
Each image has a particular geometry acquisition and does not repeat ex- 
actly due to the differences in orbital coordinates (Lauknes, 2004). The dif- 
ference phase information obtained from the sensor ASAR (in two different 
positions of the passage) is essential information to generate the i nterfero- 



gram (Lillesand et al. 2004; J arisen, 2000; Mather, 1999, Cracknel I and 
Hayes, 1991; Ferreti et al. 2,006; Farina et al. 2006). Assuming the phase 
and for a particular pixel, the acquisition identical Om and Os (Figure 1), 
the main value of the interferometric phase is expressed in Equation (1), 
where dR is the difference in slant range for a given pixel, and the length A 
wave sensor (Sandwell and Price, 1998). 




GEODETIC REFERENCE SYSTEM 1930 



Figure L Basic Geometry Dl nSAR acquisition dual pass. 



As the acquisition conditions between Om and Os , other terms that change 
the final value of the interferometric phase could be appear. The difference 
phase pixel by pixel between Pi and P2 (see Figure 1) is expressed in 
Equation (2) whereBp is the perpendicular on th e base B, A the wavelength, 
the height difference dh between Pi and P2; dp is ground deformation be- 
tween acquisitions, and the term s represents the noise surrounding the 



signal, atmospheric artifacts, radio broadcasting properties of targets, 
thermal noise of the sensor, among other sources (Massonet and Rabaute, 
1993). 



4tt£ 4tt5 4 
Atp = — dR + — dh + — + s (2) 

'hR.tanQ 'hR.srnQ A 

The subsidence in the ground is tied to the Geometry Dl nSAR acquisition 
geometry, and when the contribution of the first two terms of Equation (2) 
are removed (and assuming the fourth term negligible in the case of high 
temporal and spatial correlation between the images), the phase depends 
only on the deformation of the soil. The cyclic property of the interferomet- 
ric phase (0 to 2n) of each pixel requires it be removed to analyze the sub- 
sidence of soil, within the process called unfolding phase in calculating the 
absol ute phase of pi xel s ( M ati as, 2006) . 



3. Some studies and applications of techniques Dln- 
SAR PSInSAR in Brazil 

In this section, some studies of Brazilian experience with PSInSAR in public 
domain applied to real cases will be shown. 

In Paradella et al. (2005) present the use of SAR data and prospects for 
geol ogi cal mappi ng appl i cati ons and wi th the i mmi nence of radar operati on 
multi polarized airborne, polarimetric and interferometric SIVAM-SIPAM 
and prospects for the ALOS/ PALSAR and RADARSAT-2. 

Some multi disciplinary Brazilian works using SAR images application can 
be found in:. A) Mura (2005) and Gomes et al., (2005) have a forest appli- 
cation; B) deforestation and fire are presented in Bufalo and Valeriano 
(2005); C) Geomorphologyin France et al., (2005) and Elmiroetal. (2005) 
within the I nSAR technique. Another application involving radar images is 
shown in D) Pi ceo et al. (2007) which addresses the development of filter 
for SAR images, in Geomorphology (Dutra et al. 2005), and altimetry 
(Oliveira et al., 2009). The foil owing sums up what has been developed in 
these articles. 

I n the work developed by M ura (2005) have been two methods for i improv- 
ing the consistency with projection vectors to represent the backscatter sig- 
nal of the targets. It is used in this article polarimetric and interferometric 



data of banda P, acquired in theTapajos National Forest in 2000. It is ob- 
served that due to the penetrati on of the mi crowave si gnal i n the P band has 
not only the acquisition of interferometric phase practically at ground level, 
but also parameters that describe the structure of the forest. 

Gomes et al. (2005) show different processing digital images made to rat- 
ings in J ERS-3/SAR SPRING environment and evaluate these results from 
the Kappa coefficient, calculated from the confusion matrices. We also ana- 
lyzed the results using filters before the process of segmentation and classi- 
f i cati on of these i mages. 

Buffalo and Valeriano (2005) show in their work that can obtain verifica- 
tion of areas burned detectable from the data mosaic of J ERS-1 (SAR) of 
GRFM. They check that from the refinement of the methodology has a 
complement of PRODES Digital, developed bylNPE. 

France et al., 2005 show in his work an assessment of the loss of infor- 
mation in the mapping of fluvial forms due to the degradation of spatial 
resolution in an image RADARSAT.C band. The area chosen for the job is 
set along the lower Amazon, the southern city of Obi dos, in the state of Para 

Elmiroetal. (2005) developed his work on digital elevation models (DEM) 
generated from SAR i nterferometry obtained in the Amazon region for a 
correction method based on MDE surfaces errors determined by tri angula- 
tion with linear interpolation. The three-dimensional models in use were 
obtained from a system of SAR polarimetry complete banda P and H H po- 
larization in theX band AeroSensing RadarSysteme GmbH held that flight 
in theTapajos National Forest and vicinity. 

Picco et al. (2007) present a brief review of a multiplicative model for SAR 
images univariate discusses new filters MAP, describe the experience of 
Monte Carlo and criteria to evaluate the performance of the filters present- 
ed. 

Oliveira et al. (2009) show the initial results achieved in the generation and 
evaluation of a MDE altimetricTerraSAR-X generated from a stereo pair of 
i mages acqui red in SpotLight mode(SL). 

Guerra, Cecarelli and Lomonaco (2011) studied the viability of the Dlnsar 
and PSI nSAR with TerraSAR X to monitoring the surface around the Salto 
Pilao-SC hydroelectric. 

Finally, Paradella et al. (2012) show a review of the main SAR and I NSAR 
concepts and application perspectives of these technologies in Brazil. 



3.1. Area of interest 

The neighborhood of BoaViagem, Recife, Pernambuco, Brazil (Figure 2) is 
located south of Pernambuco Shear Zone, the main geological structure and 
structural control that divides the basin sedimentary Paraiba and Pernam- 
buco Basin. This neighborhood has overexploitation of aquifers sectors and 
to the north, a large wetland vegetation - Mangrove Park. Boa Viagem has 
variable range of sand with palm trees, liquid surfaces on the edge of the sea 
and mangroves surfaces whose dielectric parameters and polari metric sig- 
natures generate temporal decorrelation and should be excluded from pro- 
cessing images. This area is geologically sedimentary and overpumping of 
wells reduces the pressure exerted by water generating support redistribu- 
tion of grains and soil subsidence, which can reach a few centimeters per 
year and extendi ng several kilometers (Cabral et al., 2006). 

The site http:/ / www.abequa.orq.br/ trabal hos/ qerenci amento_ 363. pdf 
shows some results about aquifer monitoring of the study area. 




Figure 2. Neighborhood of Boa Viagem, Recife, Brasil. (a) Location of the area. 
Source: Santos (2005) . (b) Photo of the urban beach of BoaViagem. Fonte: Teresa 
Maia- www.pernambuco.com 



The temporal decorrelation depends on the stability of the back-scatteri ng 
of targets between the acquisiti on dates of the i mages. 

Boa Viagem in the buildings of concrete and metal fixed targets be appro- 
priate to have good correlation (good retro- reflectors in C band. Decorrela- 



tion generate targets susceptible to deletion or change of position, orienta- 
tion and shape, as trees, fences, vehicles (heavy traffic in Recife), mounts 
regular temporary stage on the beach typical of Boa Viagem, paving roads, 
and any other built structure removed or not present at all radar images. 
Other targets can generate decorrelation as construction of new buildings 
within very dynamic removal of houses and old buildings replaced by tow- 
ers of 30 floors in Recife, installation of guard rails, walkways, electronic 
speed bumps, bridges, viaducts, suspended pipelines, networks electricity 
transmission, etc.. The local knowledge of these urban characteristics is 
important in understanding the temporal decorrelation during the image 
processing and analysis and appropriate choice of spectral filters or com- 
puter algorithms to be applied. 

I n the images ASAR (Advanced Synthetic Aperture Radar) of EN VI SAT (C 
band wavelength of 3.8 cm to 7.5 cm - frequency range from 4.0 to 8.0 
GHz) modes swath IS2 and I S3 in orbits descending and ascending cover 
such additional points on the ground in the neighborhood of Boa Viagem, 
as the pattern and road corridors between buildings, which are referred to 
the most appropriate measurement of soil subsidence. The combination of 
the mentioned orbits seeks to circumvent effects of ASAR images on flat 
terrain with many buildings generating vertical radar shadows caused by 
east- west line north-south buildings very close together points covering the 
soil capable of GNSS monitoring, complicating the comparison of GNSS 
and altimetry altimetry with DlnSAR technique. 



4. Conclusion 

The overexploitation of groundwater with the use of wells for situations 
requiring a demand exceeding the avail ability of the subterranean reservoir, 
the well tends to lead to exhaustion altering the structure of the soil and 
enhancing the lowering of the soil due to compression underlying porous 
layers, a phenomenon known as soil subsidence. 

This work showed a review of work on the Dl nSAR technique with potential 
application to i nterferometry to Boa Viagem, Recife, Brazil. From this tech- 
nique it is possible to verify the changes on the ground with overexploita- 
tion. As seen DlnSAR technique is independent of other techniques for 
measuring subsidence whereas the ASAR sensor system is not located in the 
soil-building. 

Other types of phenomena (seismicity, underground works, etc.)., Besides 
the overexploitation of aquifers, subsidence can cause soil and occur in the 
same time zone and area of overexploitation of aquifers, maki ng the separa- 
ti on of causes of soi I subsi dence. 



The consideration of the geological characteristics of Boa Viagem and the 
current I and use and trends are essenti al for the processi ng of ASAR i mages 
using the DlnSAR technique, a technique that is adequate to the problem of 
measuring soil subsidence in Recife and depending on the his greatest sci- 
entific dissemination regarding technical PSInSAR. 
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